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We report an updated version of the ABKM09 NNLO PDF fit, which includes the most recent 
HERA collider data on the inclusive cross sections and an improved treatment of the heavy-quark 
contribution to deep-inelastic scattering using advantages of the running-mass definition for the 
heavy quarks. The ABMl 1 PDFs obtained from the updated fit are in a good agreement with the 
recent LHC data on the W- and Z-production within the experimental and PDF uncertainties. We 
also perform a determination of the strong coupling constant Us in a variant of the ABM 11 fit 
insensitive to the influence of the higher twist terms and find the value of a.; = 0.1 133(1 1) which 
is in good agreement with the nominal ABMl 1 one and our earlier determination. 
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Modern phenomenology of high-energy hadronic interactions in great extent is driven by the 
precise calculations, which include higher order perturbative QCD corrections and require accurate 
knowledge of the non-perturbative input coming from the parton distribution functions (PDFs). 
In turn the latter are steadily improved due to the progress in the calculations of the perturba- 
tive corrections, which are commonly used to tune the PDFs, and refined measurements of those 
processes. The ABKM09 PDFs [1] were also recently updated to the ABM 11 ones [2] taking 
advantage of more accurate deep-inelastic-scattering (DIS) data obtained by the HERA experi- 
ments [3, 4] and improved theoretical treatment of the heavy-quark electro-production. This im- 
provement is particularly important for the \ov/-x DIS data used in the analysis due to a significant 
contribution of the heavy quarks in the inclusive 
cross sections. In the ABKM09 fit it is computed 
within the 3 -flavor factorization scheme with the 
heavy quarks appearing in the final state. The 
ABMll fit is also based on the 3-flavor scheme, 
however, the MS -mass definition is now used for 
the charm and bottom quarks, in contrast to the 
pole-mass definition employed in the ABKM09 
fit. This approach provides improved conver- 
gence of the perturbative QCD calculations for 
the heavy-quark electro-production [5] which is 
crucial in view of the full NNLO conections to 
this process are still missing. The quality of the 
ABMll fit for the newly appended HERA data 
is quite good. The value of = 537 is obtained 
for the combined Hl&ZEUS sample [3] includ- 
ing the neutral-current and the charged-current 
cross section data with the total number of data 
points (NDP) equal to 486, while for the neutral- 
current data collected by the HI experiment in the 
high-inelasticity run [4] the value of x^/NDP = 
137/130. The latter data are particularly sensitive 
to the longitudinal structure function at small x, which is in turn sensitive to the small-^; gluon 
distribution. The 3-flavor scheme NNLO-predictions based on the ABMll, JR09, MSTW08, and 
NN21 PDFs are compared with the Fi measurements of Ref. [4] in Fig. |l[ While the ABMl 1 and 
JR09 analyses are based on the 3-flavor treatment of the DIS, the MSTW08 and NN21 PDFs are 
obtained with different variants of the variable-fiavor-number (VFN) scheme with the charm quarks 
appearing in the initial state. Therefore in order to provide a consistent comparison the MSTW08 
and NN21 predictions are computed employing the complementary 3-flavor PDF sets provided by 
those groups. With such a unification the spread of the predictions in Fig. [l] displays the difference 
in the corresponding PDFs, basically the small-x gluon distributions. While the PDFs obtained 
within the 3-flavor scheme provide a better description of the data, the use of a modeling based on 
the VFN scheme apparently leads to departures at low values of Q^. 

In any case the small-^; gluons by different groups can evidently consolidate using these data. 
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Figure 1: The data on Fi versus x obtained by 
the HI collaboration [4] confronted with the 3- 
flavor scheme NNLO predictions based on the 
different PDFs (solid line: this analysis, dashes: 
JR09 [6], dots: MSTW [7]). The NLO predic- 
tions based on the 3-flavor NN21 PDFs [8] are 
given for comparison (dashed dots). The value 
of for the data points and the curves in the 
plot rises with x in the range of 1 .5 45 GeV^. 
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Figure 2: The charm- (left) and the bottom-quark (right) PDFs obtained in the fit: The dotted (red) Hnes 
denote the ±1(7 band of relative uncertainties (in percent) and the solid (red) line indicates the central 
prediction resulting from the fit with the running masses taken at the PDG values. For comparison the 
shaded (grey) area represents the results of ABKM09 [1]. 



The DIS HERA data are sensitive to the value of the charm-quark mass rUc. Its MS-value is ob- 
tained by mdrnc) = 1.24 ±0.03(exp)+o;o^(scale)+o;Q7(theory) GeV [9] preferred by the inclusive 
HERA data of Refs. [3,4] combined with the semi-inclusive data on open charm production [11] 
is in good agreement with the world average [10]. This provides an additional argument in fa- 
vor of the 3-flavor description of DIS. Furthermore, in the ABM 11 analysis the accuracy of the 
small-;c 4- and 5-flavor heavy-quark PDFs is improved as compared to the ABKM09 one due to 
the constraints on charm- and bottom-quark masses because the world average can be imposed, cf. 
Fig. ^ In contrast, the values of extracted from the HERA data within different variants of the 
VFN scheme demonstrate a wide spread and, moreover, they are systematically lower than the ones 
obtained from the e^e^ data [11]. 

The value of the strong coupling constant a^{Mz) =0.1134(8) obtained intheNNLO ABMll 
fit is in agreement with the earlier one of ABKM09. Meanwhile the NMC and SLAC DIS data sets 
prefer smaller and bigger value of a^, respectively, cf. Fig. ^ Note that these two data sets are 
sensitive to the contribution of the higher twist terms [2, 12], while the HERA and BCDMS data do 
not suffer from this contribution due to larger typical values of the momentum transfer involved. 
To avoid the influence of higher twist terms we perform a variant of the ABMl 1 fit with the SLAC 
and NMC data dropped and the higher twist terms set to zero. The value of as{Mz) = 0.1 133(11) 
obtained in this way is in a good agreement with the nominal value of ABMl 1 that demonstrates a 
consistent treatment of the higher twist terms in the present analysis. The NNLO ABMl 1 value of 
OCs(Mz) is essentially lower than the PDG world average [10]. On the other hand the ABMl 1 value 
of as(Mz) is in agreement with the recent NNLO determinations obtained from the analysis of 
thrust distributions measured in the e^e^ collisions with account of the power corrections [13, 14] 
and also the lattice calculation of the QCD static energy [15]. 

Benchmarks for the rates of key-processes at Tevatron and the LHC were given in [1,2, 16]. 
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Figure 3: The j^-profile versus the value of as{Mz), for the separate data subsets, all obtained in variants of 
the present analysis with the value of fixed and all other parameters fitted (solid lines: NNLO fit, dashes: 
NLO fit). 



For the first LHC data on the Drell-Yan differential cross sections [17-20] we observe good agree- 
ment with the predictions based on the ABM 11 PDFs. These data were obtained by detecting 
the leptonic decays of W- and Z-bosons with the charged lepton pseudo-rapidities in the range 
of -2.5-^2.5, -2.4-^2.4, and 2-^4.5 for the ATLAS, CMS, and LHCb experiments, respec- 
tively. The data on the Z-boson pseudo-rapidities and on the charged lepton pseudo-rapidities for 
the case of W-boson production are confronted to the NNLO predictions obtained with ABM 11 
PDFs in Figs. The central values of the predictions are obtained using a fully exclusive code 
DYNNLO 1.3 [21]. However, the uncertainties due to PDFs are computed with another fully ex- 
clusive code FEWZ 3.1 [22], which allows to estimate the variation of the cross section with the 
PDFs without accumulating the Monte-Carlo integration errors in the cross section differences. 
The typical accuracy of the DYNNLO results displayed in Figs. Q-^ is better than 1% and it does 
not exceed the data errors. 

Agreement between the data and the predictions is quantified by the following functional 

x' = ny'-f^)[c-'Uyj-fi\ (1) 

with the covariance matrix 

28 

Cu = C:f + Zifi'^-fhiff-ff)- (2) 
k=l 

The first term in Eq. corresponds to the covariance matrix describing the errors in the exper- 
imental data, while the second term appears due to the PDF uncertainties. The latter contribute 
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ATLAS (7 TeV, 35 1/pb) 




Figure 4: The ATLAS data of Ref. [17] on the rapidity distribution of the charged leptons (left panel), 
(central panel), and Z-boson (right panel) in comparison with the NNLO predictions computed with the 
ABMll PDFs taking into account the uncertainties due to PDFs (shaded area). The dashed curves display 
the ABMl 1 predictions obtained with the value of as{Mz) = 0.118. 



to the covariance matrix by virtue of the differences between predictions obtained for the central 
ABMl 1 PDF set, ff\ and the ones for each of the 28 ABMl 1 PDF sets representing the PDF un- 
certainties, fj'^K The experimental covariance matrix for the ATLAS data of Ref. [17] is computed 
by 

c-''=d^joi+fj^ffj:^, (3) 

k=\ 

where a, are the statistical errors in the data combined in quadrature with the uncorrelated errors, 
s'l are the correlated systematic uncertainties representing 31 independent sources including the 
normalization, and 5ij is the Kronecker symbol. In view of the small background for the W- 
and Z-production signal all systematic errors are considered as multiplicative. Therefore, they are 
weighted with the theoretical predictions . The experimental covariance matrices for the CMS 
and LHCb data of Refs. [18-20] are available in publications and they are directly employed in 
Eq. (^ after re-weighting by the theoretical predictions similarly to Eq. (^ with the normalization 
uncertainty taken into account in the same way as for the ATLAS data. 

The values of X'^ for different LHC data samples are summarized in Table 1. Due to an incon- 
sistency in the covariance matrix for the LHCb Z-boson data [23] quoted in Ref. [19] the Z-boson 
data are not employed in the calculation of for this LHCb sample. In general the values of 
X^/NDP in Table 1 are comparable with 1 within possible statistical fluctuations. The ATLAS 
data [17] systematically undershoot the predictions, although the statistical discrepancy of the dis- 
agreement is not significant in view of still relatively large experimental uncertainties. Moreover, 
it is worth noting that the recent high-luminosity LHCb data [20] on the Z-boson production are 
in much better agreement with the ABM 1 1 predictions. In particular this is also observed for the 
rapidity range overlapping with the corresponding ATLAS data that may indicate a potential in- 
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CMS (7 TeV, 840 1/pb) 
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Figure 5: The same as Fig. ^ for the CMS data of Ref. [ 1 8] on the electron asymmetry rapidity distribution. 
The curve corresponding to the ABMll predictions obtained with the value of as{Mz) = 0.118 coincides 
with the nominal one. 

consistency between the ATLAS and LHCb normalization. The ABMll predictions demonstrate a 
certain tension with the earlier LHCb data [19] on the VK+-boson production at moderate rapidity, 
however, the overall agreement is satisfactory. If, nonetheless, the discrepancies with the ATLAS 
and LHCb data being observed are confirmed with better statistical significance, the ABMl 1 quark 
distributions at .x ~ 0.01 would have to be lowered in order to accommodate those data in a fit. 
Meanwhile, the charged-lepton asymmetry which is less sensitive to the experimental uncertainties 
is in a good agreement with the ABM 1 1 predictions both for the ATLAS and CMS data. 



Experiment 


ATLAS [17] 


CMS [18] 


LHCb [19] 


LHCb [20] 


Final states 


W+ ~>l+V 

w- -^rv 




-> jU^v 


Z e+e- 


Luminosity (1/pb) 


35 


840 


37 


940 


NDP 


30 


11 


10 


9 




35.7(7.7) 


10.6(4.7) 


13.1(4.5) 


11.3(4.2) 



Table 1: The value of obtained for different samples of the Drell-Yan LHC data with the NNLO ABMl 1 
PDFs. The figures in parenthesis give one standard deviation of equal to \/2NDP. 
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The LHC Drell-Yan data were confronted to the predictions based on the different PDFs in- 
cluding ABM 11 ones in Ref. [24] earlier. In contrast with our analysis the NNLO predictions of 
Ref. [24] are obtained using the NLO MCFM code [25] supplemented by the NNLO K-factors 
which were computed separately with the DYNNLO code. Furthermore, the benchmarking of 
Ref. [24] is performed using the PDFs corresponding to the value of a^(A/z) = 0.118. This value 
is disfavored by the ABM 11 fit. Nonetheless, in order to allow a consistent comparison with the 
benchmarking of Ref. [24], we also provide the ABM 11 prediction corresponding to the value of 
oCsi^z) = 0.1 18. The predictions obtained in this way are compared with the data and the nominal 
ABM 11 predictions in Figs. ^[,^,0. The charged-lepton asymmetry does not depend on the 
value, while the shift in the W and Z absolute cross-sections obtained with this variation of is of 
the order of the uncertainties due to PDFs. The value of = 14.3 obtained for the LHCb data [19] 
with the enhanced value of a, is within the statistically allowed range. Meanwhile the value of X'^ 
obtained in this case for the ATLAS data [17] is 41.5, i.e. slightly above the la statistical margin. 
This is obviously related to the offset of the nominal ABM 11 predictions w.r.t. the ATLAS data. 
Therefore the final conclusion about the significance of this excess can be drawn only after the nor- 
malization discrepancy for ATLAS and LHCb is clarified. Anyway, both values are significantly 
smaller than those quoted for the ABM 11 benchmarking in Ref. [24]. In particular this happens 
due to the PDF uncertainties were not taken into account in the statistical analysis of Ref. [24]. The 
central values for the ABM 11 predictions obtained in that benchmarking also do not agree with 
ours in places, cf. Fig. 12 in Ref. [24] and Fig. ^ in the present paper The recent LHCb data on 
the Z — )• e+e^ production were not included into the benchmarking of Ref. [24], however, they are 
compared with the NNLO predictions based on the CTIO [26], MSTW09 [7], and NN21 [27] PDFs 
in Fig. 4 of Ref. [20]. These PDFs provide consistent results. Nonetheless all three undershoot the 

'To facilitate further clarification of the discrepancy the tabulated NNLO ABM 11 predictions for the available 
ATLAS, CMS, and LHCb Drell-Yan data are attached to the arXiv version of these Proceedings. 
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LHCb (7 TeV, 940 1/pb) 



NNLO ABMll 
oc^(M2)=0.1134 

NNLO ABMll 
oc/M^)=0.118 




Figure 7: The same as Fig. for the LHCb data of Ref. [20] on the Z-boson rapidity distribution. 



LHCb Z-boson data at large rapidity. To quantify this disagreement we check for these data 
taking NNLO MSTW08 predictions with account of the PDF uncertainties and find the value of 
26.2, more than twice bigger than ours, cf. Table 1. 

In summary, we present the updated version of the NNLO PDF fit based on the combination 
of the DIS and Drell-Yan data. The update includes the most recent HERA collider data on the 
inclusive neutral-current and charged-current cross sections. From the theory side we improved the 
treatment of the heavy-quark contribution to DIS using advantages of the running-mass definition. 
The ABMll PDFs obtained from the updated fit are in a good agreement with the recent LHC 
data on the W- and Z-production within the experimental and PDF uncertainties. Our results do 
not confirm the benchmarking of the ABMll PDFs published elsewhere [24]. We also perform a 
determination of the value in a variant of the ABMl 1 fit insensitive to the higher twist terms and 
find a value, which is in very good agreement with the nominal one and our earlier determination 
of ttv 

Acknowledgments. We acknowledge fruitful discussions concerning the Drell-Yan LHC data 
with Uta Klein and Katharina Miiller. We are also grateful to Juan Rojo for clarification of some 
details of Ref. [24]. 
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